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ABSTRACT 

We present detailed parameter determinations of two chemically normal late A-type stars, 
HD32115 and HD 37594, to uncover the reasons behind large discrepancies between two 
previous analyses of these stars performed with a semi-automatic procedure and a "classical" 
analysis. Our study is based on high resolution, high signal-to-noise spectra obtained at the 
McDonald Observatory. Our method is based on the simultaneous use of all available observ- 
ables: multicolor photometry, pressure-sensitive magnesium lines, metallic lines and Balmer 
line profiles. Our final set of fundamental parameters fits, within the error bars, all available 
observables. It differs from the published results obtained with a semi-automatic procedure. A 
direct comparison between our new observational material and the spectra previously used by 
other authors shows that the quality of the data is not the origin of the discrepancies. As the 
two stars require a substantial macroturbulence velocity to fit the line profiles, we concluded 
that neglecting this additional broadening in the semi-automatic analysis is one origin of dis- 
crepancy. The use of Fe i excitation equilibrium and of the Fe ionisation equilibrium, to derive 
effective temperature and surface gravity, respectively, neglecting all other indicators leads to 
a systematically erroneously high T e ff. We deduce that the results obtained using only one 
parameter indicator might be biased and that those results need to be cautiously taken when 
performing further detailed analyses, such as modelling of the asteroseismic frequencies or 
characterising transiting exoplanets. 
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1 INTRODUCTION 

The advent of space missions aiming to obtain very accurate pho- 
tometry for an increasing number of stars (e.g. CoRoT and Ke- 
pler) led to the necessity of a large scale work to obtain high pre- 
cision spectroscopic fundamental parameters, effective temperature 
in particular, to allow i.e. the modelling of the pulsation frequencies 
or the characterisation of transiting exoplanets. This large spec- 
troscopic analysis campaign can be performed within reasonable 
time-scales only with automatic and semi-automatic procedures. It 
is therefore crucial to critically compare some of the results ob- 
tained in this way with other independent, more "classical", meth- 

* Data obtained with the 2.7-m telescope at McDonald Observatory, Texas, 
US. 
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ods and highlight the discrepancies to understand their origin and 
imp rove these procedure s. 

iBruntt et alj d2010l) derived fundamental parameters for a set 
of 23 solar-type stars adopting various different techniques (e.g. 
asteroseismology, parallax, spectroscopy), concluding that purely 
spectroscopic methods lead to results comparable to the more ro- 
bust ones, although small corrections might be necessary. It is im- 
portant to carefully assess whether there are cases where the spec- 
troscopy fails in recovering the correct set of parameters, and if 
such erroneous parameters are produced, we need to understand 
why and correct the methodology. 

In this work we concentrate on a few discordant results ob- 
tained with the semi-autom atic procedure developed by H. Bruntt 
(see e.g. lBruntt et al J2002h and adopted by many authors to analyse 
several different types of stars. 

For the solar-type pulsator HD 49933, iGillon & M again 
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d2006h and lBruntt etall j2008l) (hereafter B08), making use of two 
very similar semi-automatic proc edures, derived an effectiv e tem- 
perature (r e jf) of about 6750 K. iRvabchikova et all l l2009L here- 
after R09) re-analysed this same star applying a "classical" anal- 
ysis (based on photometry, equivalent widths and hydrogen lines) 
to higher quality spectra, obtaining a r e j_of_6500K, demonstrat- 
ing this is a better value for T sff . Later iBrunt 3 fe009h reanalysed 
the same spectra of HD 49933 used by R09 with the same method 
used by B08 obtaining a r eff of 6570 K, in agreement with R09, 
but in disagreement with B08. The later analysis did not explain 
the reason for such disparate results from the same semi-automatic 
procedure on the same star, leaving open the possibility that differ- 
ent results might be obtained with data from different instruments. 
This possibility needs to be ruled out and it is important to un- 
derstand the origin of these discrepancies, furthermore because the 
code employed by B08 and iBrunttl f2009h is adopted for most of 
the spectroscopic analysis of the C0R0T stars. 

To understand the discrepancies, we must examine more than 
one star, therefore we decided to extend our critical analysis of 
B08's methodology for fundamental parameter determination to 
two oth er stars: HD 321 15 an d HD 37594, used as comparison stars 
by B08. Bikmaev et al. (2002) showed the results of a fundamental 
parameter determination and abundance analysis of these two stars 
performed with a classical method, obtaining r c(f of 7250 K and 
7170 K, respectively. B08 re-analysed HD 32115 and HD 37594 
employing a semi-automatic procedure on a different set of spec- 
tra, obtaining 7670 K and 7380 K, respe ctively. Here we use h igher 
quality data than that previously used by Bikma ev et al. IHH) and 
B08, focusing on the parameter determination and using only care- 
fully selected spectral lines considering their log gf values, blend- 
ing, and known non-LTE effects affecting those lines. 

HD32115 is a single line spectroscopic binary for which 
iFekel et alj (2006) determined an orbital period of about 8 days, 
and concluded that the companion is either a late K- or an early 
M-type star. This ensures that the spectral lines of the primary are 
not affected by the companion at any wavelength and that it is also 
safe to perform a spectrophotometric analysis to derive the stellar 
parameters of the pri mary star (s e e Sect . l4.1.2t . 

For HD 32115, TFekeletalJ {2006) derived an effective tem- 
perature of 725 1 K and a surface gravity of 4. 26. From the HIP - 
PARCOS parallax and from their orbit analysis, Fekel etal.U2006l) 
derived a stellar mass and radius of 1.5 M/M e and 1.5+0.1 R/R a , 
respe ctively. By means of stellar struc ture and evolution calcula- 
tions, lAUende Prieto & Lamber 3 dl999h derived an effective tem- 
perature of 7413 K and a surface gravity of 4.29+0.04. They de- 
rived also a stellar mass of 1. 54± 0.06 M/M n and a st ellar radius 
of 1.48±0.03^/i? Q . We notice that lBikmaev et al.1 J2002h 's T cff for 
HD32115 c oinci des with the effective temperatures obtained by 
IFekel et alj J2OO6I) . On the other hand, B08's r eff does not agree 
with any of the previous determinations. 



2 OBSERVATIONS 

Spectra of HD32115 and of HD 37594 were obtained on 2010, 
November 30th with the Robert G. Tull Coude Spectrograph (TS) 
attached at the 2.7-m telescope of McDonald Observatory. This is 
a cross-dispersed echelle spectrograph yielding a resolving power 
(R) of 60 000 for the configuration used here. The signal-to-noise 
ratio (S/N) per pixel at A -5000 A is 490 and 535, respectively for 
HD32115 and HD 37594. 

Bias and flat-field frames were obtained at the beginning of 



the night, and a Th-Ar spectrum, for wavelength calibration, was 
obtained between the two science spectra, which were re duced 
using the Image Reduction and Analysis Facility (IRAlfl, iTodvl 
1993). Each spectrum, normalised by fitting a low order poly- 
nome to carefully selected continuum points, covers the wavelength 
range 3633-10849 A, with several gaps between the orders at wave- 
lengths greater than 5880 A. Making use of the Th-Ar spectrum we 
ensured the stability of the resolution. The CCD count rates were 
well below the saturation level, ensuring linearity and therefore no 
systematic differences between strong and weak lines. This was 
also confirmed by the complete absence of any correlation between 
line abundance and wavelength (see Sect |3.2.31 . 

The normalisation of the Hy line was of crucial importance 
since we adopted the profile fitting of the Hy line wings as one of 
the temperature indicators. We were unable to use either Ho- and 
H/3 because the former was not covered by our spectra and the 
latter was affected by a spectrograph defect, preventing the nor- 
malisation. We were able to perform a reliable normalisation of 
the Hy line us i ng the artificial flat-fielding technique described by 
Barkl enTet al.1 j2002h . This normalisation procedure was already 
successfully applied t o data obtained with this spectrograph by 
iKolenberg et al.1 fcOld) . 



2.1 Comparing spectrographs 

As mentioned in Sect. □ Bruntt(2009) left open the possibility that 
the discrepancies obtained for HD 49933 with their previous re- 
sults (B 08) were caus ed by systematic differences in the observed 
spectra terunttll200^ . analysed HARPS spectra, while B08 anal- 
ysed CORALIE spectra, where both instruments operate at LaSilla, 
but on different telescopes). We checked if this is the case. This 
control is important, to remove the observed spectra as a possible 
source of systematic uncertainties and also to check the quality of 
the normalisation, as independently performed on spectra obtained 
with different instruments. 

Here we compare the spectra obtained with TS and CORALIE 
(used by B08), which is an echelle spectrograph, reaching a res- 
olution of R ~50000, mounted on the 1.2-m Euler telescope in 
La Silla, Chile. Details of the CORALIE spectrog raph and on the 
data reduction can be found in lDe Cat et al.1 (2006). 

Figure [JJ shows a comparison between the TS and the 
CORALIE spectrum of HD 321 15 in the wavelength range around 
the strong Fen /15018 A line. In this plot we display also the dif- 
ference spectrum (in % shifted upwards of 0.9) between the two 
spectra. 

The rms of the difference spectrum (a portion is shown in 
Fig. [TJ, calculated on several continuum regions is comparable to 
the rms of the CORALIE spectrum. We obtained the same conclu- 
si on comparing the TS and CORALIE spectra with the ones used 
bv lBikmaev et alj J2002h . The spectra of HD 37594 demonstrate an 
identical behavior. 

These comparisons let us conclude that there are no significant 
differences b etween the spectra us ed in the present work and those 
employed by Bi kmaev et alj 12002J) and by B08, thus excluding the 
quality of the data as the origin of the discrepancies described in 



1 IRAF (http://iraf.noao.edu/) is distributed by the National Opti- 
cal Astronomy Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation. 
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Figure 1. Comparison between the spectra of HD 32115 obtained with the 
CORALIE spectrograph (dashed red line) and the TS spectrograph (thick 
black line) in the wavelength region around the /15018 A Fen line. The dif- 
ference (in %) between the two spectra is shown by the dashed- dotted blue 
line. The difference spectrum is shifted upwards by 0.9. No difference is 
visible between the two spectra, except for the core of the strong Fe n line 
at /15018 A which displays the difference in resolution between the two in- 
struments. 



Sect.Q] In addition, this comparison sets an upper limit of ~1% on 
the uncertainty due to the normalisation. 



3 THE FUNDAMENTAL PARAMETERS 

We adopted photometric indicators to set the starting point in the 
determination of the fundamental parameters, which we refined 
making use of hydrogen lines, metal lines, and, as a final check, 
synthetic colors and the spectral energy distribution. Spectroscopic 
tuning of the fundamental parameters is needed because different 
photometric systems and calibrations give different parameters and 
uncertainties. Spectroscopic analysis, performed with in this way, 
will produce a parameter set which fits all the indicators consis- 
tently, the uncertainties are thus reduced for those derived from 
photometric analysis alone. 

We computed model atmospheres of HD 32115 and HD 37594 
using the LLmodels stellar model atmosphere code dShulvak et al.l 
2004). For all the calculations Local Thermodynamical Equi- 
librium (LTE) and plan e-parallel geometry were assumed. We 
used the VALD database dPiskunov et alJll995l:lKupka et alJI 19991 ; 
iRvabchikova et all 1 19991) as a source of atomic line parameters 
for op acity calculations. Convect i on wa s implemented according 
to the Canuto & Mazzitelli ( 1991, 1992) model of convection (see 
iHeiter et alfeOOl for more details). 



3.1 Photometric indicators 

Initial guesses for the effective temperature (r c ff ) and surface grav- 
ity (log g) were obtained from calibrations of different photometric 
indices for normal stars. The effective t emperature and gravity were 
derived from Stromgr en photometry i Hauck & M ermilliod 19981) 
with c a libration s by Moon & Dworetskv 1 1985h, Napiwotz ki et al.l 
dl993h . iBaional dl994l) . iRibas et al.l dl997h. ICastelli et alj dl997l) . 
and from Geneva photo metry dRufenedl 1988b with calibrations by 
lNorth&NicoTetldl990l) . 

Table[T]summarises the set of r cff and log g obtained with each 
adopted calibration and photometry, reflecting the scattering due to 
the use of different calibrations and photometric systems. 



The interstellar reddening plays an important role when con- 
verting the observed photometry into fundamental parameters. We 
calculated the interstel lar reddening adopting th e galactic extinc- 
tion maps provided bv lAmores & Lepi ne (2005), obtaining E(B - 
V)=0.00 for both stars, which we adopt for the determination of the 
synthetic colors (see Sect. 14.1. It and spectral energy distribution 
(see Sect. 14.1.21 . This is also in agreement with ot her models of in- 
terstel lar extinction in the solar neighborhood, e.g. lLallement et all 
d2003h . 

Excluding the results obtained with the calibration by 



ICastelliet all d 19971) . which gives much lower T e ff and log g com- 
pared to the others, we set the center for the calculation of our 
model grid to T cff = 7250 K and logg = 4.2, for HD32115 and 
to T eff = 7100 K and logg = 4.2, for HD 37594, adopting steps of 
50 K in r eff and 0. 1 in log g. 



3.2 Spectroscopic indicators 

3.2.1 Hydrogen lines 

For a fully consistent abundance analysis, the photometric param- 
eters must be checked and eventually tuned according to spectro- 
scopic indicators, such as hydrogen line profiles. In the temperature 
range where HD 321 15 and HD 37594 lie, the hydrogen line wings 
are sensitive almost exclusively to T e g variations. To spectroscopi- 
cally derive r eff from hydrogen lines, we fitted sy nthetic line pro- 
files, calculated with SYNTH3 dKochukhovfeoOTI) . t o the observed 
Hy pr ofiles. SYNTH3 incorporates the code by iBarklem et al.l 
(20000 that takes into account not only self-broadening but also 
Stark broadening (see their Sect. 3). For the latter, the default mode 
of SYNTH3, ado pted in this wo rk, uses an improved and extended 
HLINOP routine dKuruczlll993h . 

Figure [2] shows a comparison of the observed Hy line profiles 
of HD32115 and HD 37594, with two synthetic profiles for each 
star. One synthetic profile was calculated with our final set of pa- 
rameters while the other one with the set of parameters by B08. 

For both stars, the synthetic spectrum corresponding to our fi- 
nal model (see Sect. |4) does not fit the Hy profile perfectly, this 
would require a lower temperature by -80-100 K. We attribute 
most of the difference between our final synthetic and observed 
Hy profiles to the normalisation, which is always challenging for 
hydrogen lines observed with echelle spectra. To quantify the un- 
certainty introduced by the normalisation, we compared the TS Hy 
profile with the one we obtained from the CORALIE spectrum. The 
maximum difference between the two profiles (independently nor- 
malised) is ~2%, less than the difference introduced in the synthetic 
spectrum by changing T e g by 100 K. 

Figure|2]shows that the effective temperatures adopted by B08 
for the two stars are too high to even remotely fit the hydrogen line 
profiles. 



3.2.2 Gravity from metallic lines with extended wings 

The surface gravity was derived from two independent methods 
based on line profile fitting of Mgi lines with developed wings 
(analysis described here) and ionisation balance for several ele- 
ments (an alysis performed in Sec t. I3.2.3t . The first method is de- 
scribed in lFuhrmann et al.l <jj997) and is based on the fact that the 
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Table 1. Set o f T e ff and logg obtained for HD32115 and HD 37594 using different calibrations for Stromgren <Hauck & Mermilliod[T993) and Geneva 
lRufenedfl988l) photometry. The las t two lines give the average values of T e g and log g, with their uncertainties, calculated without taking into account the 
results from the Castelli et al. ( 1997) calibration. 



HD32115 HD 37594 



Photometry 


r eff [K] 




r eff [K] 


log? 


Calibration 


Stromgren 


7296 


4.23 


7171 


4.22 


Moon & Dworetskv (1985) 




7207 


4.13 


6899 


3.86 


Napiwotzki et al. (1993) 




7421 


4.11 


7276 


4.05 


Balona (1994) 




7308 


4.27 


7189 


4.32 


Ribas et al. (1997) 




6998 


3.84 


6667 


3.47 


Castelli et al. (1997) 


Geneva 


7263 


4.47 


7157 


4.47 


North &Nicolet (1990) 




7300±80 


4.24±0.14 


7140±142 


4.18±0.24 
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Figure 2. Comparison between the observed Hy line profile (black solid 
line) and synthetic profiles calculated with our final adopted parameters 
(red dashed line) and with the parameters adopted by B08 (blue dotted 
line). The upper comparison is for HD 37594, while the lower comparison 
is for HD 32115. Our final adopted set of parameters is T^g = 7250± 1 00 K, 
logg = 4.2±0.1 for HD 32115 and T e g = 7150±100K, logg = 4.2±0.1 for 
HD 37594. B08's set of parameters is 7/ cff = 7670±170 K, logg = 4.44±0.13 
for HD 321 15 and r cff = 7380± 190 K, log g = 4.08±0. 16 for HD 37594. The 
profiles shown for HD 37594 were rigidly shifted upwards by 0.5. 



wings of the Mgi lines at AA5167, 5172 and 5183 A are very sen- 
sitive to log g variations. In practice, first the Mg abundance is de- 
termined from other Mg i lines without developed wings, such as 
A 5528 A, and then the wings of the three lines listed above are 
tuned to the log g value. 

It was not possible to use the Mg i lines at AA 5167 and 5183 A 
because the former is blended by a strong Fe i line, while the latter 
falls at the border of two overlapping spectral orders and the nor- 
malisation was not good enough for the required precise fitting of 
the line wings. In the end, we based our fit for the log g determina- 
tion just on the Mg i line at A 5 172 A. 

For HD 32115, we used the Mg i line at A 5528 A and the Mg n 
line at A 7877 A to measure the Mg abundance, while for HD 37594 
we used only the former because the latter is blended with a telluric 
line, making the abundance determination unreliable. 

Given the presence of non-LTE effects for Mg , in particular for 
the Mgn infrared line dAbia & Mashonkinal [2004 ) , we calculated 



the non-LTE corrections for the Mg lines used. As our synthetic line 
profiles are in LTE, non-LTE calculations are necessary to make 
sure that the Mg abundance, we derive from the Mg i A 5528 A and 
Mg ii A 7877 A lines, is correct and therefore applicable for the line 
profile fitting of the Mg i A 5 172 A line. 

Non-LTE corrections for neutral and singly-ionised mag- 
nesium were carried out usi ng the codes D ETA IL and SUR- 
FACE, originally developed bv lGiddingsl dl98lh and lButled d 19841) 
along with the model atmosphere computed with LLmodels. Our 
calculations take into account the recent improvements in the 
atomic data for Mg, the extensive descript ion of the model atom , 
and non-LTE line formation presented by Prz vbilla et alj d200ll) . 
For both stars the non-LTE abundance correction A(Mg/N wl ) = 
log(Mg/N m ) LTE - log(Mg/N m ) NLTE is +0.02 dex for the Mg i lines 
at ^5172 and 5528 A, while for the Mg n infrared line at A 7877 A 
w e obtained a correction of -0 .07 d ex, in agreemen t with t he results 
bv lAbia & Mashonkinal d2004l) and lPrzvbiila et alj d200ll) . 

In HD32115 for the Mgi /15528A line we obtained 
a LTE abundance of log(Mg/N, ol ) LTE = -4.50, therefore 
\og(Mg/N, ol f LTE = -4.48. Similarly, for the Mg n line at A 7877 A 
we obtained a LTE abundance of log(Mg/N lol ) LTE = -4.39, there- 
fore log(Mg/N m ) NLTE = -4.46. Spectral synthesis in the region 
of the weaker Mg n lines at AA 4390 and 4427 A, are not sensitive 
to non-LTE effects, requires an abundance in close agreement with 
that derived from the Mg i lines. Therefore to perform the fit of the 
line wings of the Mg 1/15172 A line, we set the Mg abundance for 
HD 321 15 at -4.48, in non-LTE, and -4.50, in LTE. 

For HD 37954, non-LTE corrections are similar to those in 
HD32115, in particular they are identical for both Mgi lines at 
A 5528 A and A 5172 A. Therefore we applied for the fitting of the 
Mg i A 5 172 A line wings the LTE Mg abundance derived from the 
A 5528 Aline: log(Mg//V„„) ir£ = -4.77. 

To derive logg from the fit of the Mgi lines with extended 
wings, very accurate logg/ values and Van der Waals (logy Waals ) 
damping constants are needed. Two sets of log gf laboratory mea- 
surements for the Mg i triplet are avail able. The first one ba sed on 
the lifetime laboratory measurements dAnderson et alJ[T967l) came 
from the VALD database, while the second set, based on life- 
times and branching r atio measurements was recently published by 
Alde nius et "ai] d2007h . The accuracy of this set of transition proba- 
bilities is <r(loggf)= +0.0 4 dex. Van der Waals damp ing constants 
in VALD are calculated by Barkl em & O'MaraldTo OO). Another es- 
timate of (log ywaais) was given bv lFuhrmann et alj dl9971) . who de- 
rived logywaais - -7-4 2 from the fitting of the solar lines us i ng the 
lAnderson et alj dl967l) oscillator strengths. Fuhrmann et alj dl997t) 
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noted that Stark broadening does not practically influence the Mg i 
line profiles in the solar spectrum and in Procyon, while it might 
be more significant in hotter stars. For our analysis we employed 
Stark damping constant logy Slark = -5.44 for the Mgi triplet 
and logystark = -4.63 for the Mg i A 5528 A line calculated by 
iDimitriievic & Sahal-Brechou dl994l) . Our calculations show that 
the line profile of the latter line in the spectra of both HD32115 
and HD 37594 is not sensitive to Stark and Van der Waals broaden- 
ing effects. 

First, we checked the ato mic parameters on the NSO solar 
flux atlas dKurucz et alj Il984h . LTE synthetic spectrum calcula- 
tions in the region of the Mg i triplet and of the Mg i A 5528 A 
line were perform ed for three different m odels of the solar atmo- 
spher e: MARCS dGustafsson et alj|2008l) . MAFAGS dGrupp et all 
2009), and the one calculated with the LLmodels code. We fit the 
extended wings, not the cores of these lines which are sub|ect to 



non-LTE effects. Using logywaais from Barklem & O'Mara 
for all lines and transition probabilities from lAldenius et al 



(2000) 



(2007) 



we derived the following LTE Mg abundance in the solar atmo- 
sphere: log(Mg/7V,ot) = -4.54 (MARCS), -4.55 (MAFAGS), and 
-4.51 (LLmodels). It corresponds to 7.50, 7.49 and 7.53 in logarith- 
mic scale where \og(H)= 12.00. Our estimates are closer t o the Mg 
meteoritic abundance of 7.53±0.01 dLodders et aT] |2009). then to 
the most recent publishe d value of the M g abundance of 7.60+0.04 
in the solar photosphere ( Asplund et al. 20091). 

For HD32115 and HD 37594, careful fit of the Mgi 
/i 5 172 A lin e profile, calc ul ated with the transition prob- 
abili ties by | Aldenius et"al] d200"7|) and damping constants 
fromlBarklem & O'Marij [2000 ) and lDimitriievic & Sahal-BrechoJ 
dl994n T results in the final value of logg=4.2±0.1, in good agree- 
ment also with the gravity estimates from the photometric cali- 
brations. Error estimates include the claimed ±30 % error in the 
Stark damping constant calculations and also the possible errors in 
other spectral line parameters. The log g value of 4.2±0 .1 is also in 
very g ood agreement with the results obtain ed by both lFekel et al.l 
d2006l) and lAllende Prieto & Lambertl d 19991) for HD 321 15. 

Figuref3]shows comparisons for HD 321 15 and HD 37594 be- 
tween the observed and synthetic line profile of the A 5 172 A Mgi 



line, c alculated adopting the oscill ator strength fro m Aldeni us et al.l 
( 2007) and damping constant from lBarklem & O'Marald200d) . The 
synthetic profiles calculated with the parameters and abundances by 
B08 are also shown (B08 adopted oscillator strength and damping 
constant given in VALD). 

Figure [3] shows an excellent agreement between the synthetic 
profiles calculated with our final adopted parameters and the ob- 
served spectrum. On the other hand, the synthetic profiles calcu- 
lated with B08's parameters display a disagreement with the obser- 
vations, in particular for HD 321 15. For HD 37594 the discrepancy 
is rather small: the effect of the higher is mostly compensated 
by the lower gravity. 

The other two Mgi lines of the triplet (,U5167 and 5183 A) 
will provide the same results we obtained with the A 5 172 A Mgi 
line, as tested by R09 with the Sun, Procyon and HD 49933. 

For HD 321 15, we use d the stellar mass and bolometric mag- 
nitude bv lFekel et alj 12006J), in addition to our r cff , to derive log g: 
log(g/go) = log(M/M G ) + 4 • log(7W7V /)0 ) + 0.4 • (M bol - 4.75). 
We obtained a logg value of 4.24, in goo d agreement with our 
previous estimation and w i th the results bv lFekel et alj d2006h and 
lAUende Prieto & Lambertl d 19991) . 




Wavelength (A) 



Figure 3. Comparison between the observed Mg i line profile of the line at 
A 5172 A (black solid line) and synthetic profiles calculated with our final 
adopted parameters (red dashed line) and with the parameters adopted by 
B08 (blue dotted line). The upper comparison is for HD 37594, while the 
lower comparison is for HD 321 15. The profiles shown for HD 37594 were 
rigidly shifted upwards of 0.3. Our synthetic profiles a re calculated adopt- 
ing t he oscillator strength from lAldenius et alj |2()()7|. damping constant 
from lBarklem & O' Mara (2000) and a LTE Mg abundance of -4.50 dex for 
HD 321 15 and -4.77 dex for HD 37594. B08 does not provide a Mg abun- 
dance for HD 37594, therefore we derived it from the Mg i line at A 5528 A, 
adopting the parameters given in their paper. 



3.2.3 Metallic lines 

The metallic-line spectrum provides further constraints on the at- 
mospheric parameters. If no deviation from LTE is expected, there 
should be no trend in individual line abundances as a function of 
excitation potential. Examining this for any element/ion therefore 
provides a check on the determined value of r cff . The balance be- 
tween different ionisation stages of the same element similarly pro- 
vides a check on logg. The microturbulent velocity (v mlc ) is sub- 
sequently determined by minimising any trend between individual 
abundances and equivalent widths for a certain ion. Determining 
the fundamental parameters in this way must be done iteratively 
since, for example, a variation in r c( j leads to a change in the best 
logg and u mic . This methodology of fundamental parameter deter- 
mination from the metallic line spectrum is adopted i n almost all 
semi-automatic abundance analysis procedures (e.g. iBruntt et al.l 
l2002l : ISantos et alj2004lGillon & Magair]|2006l) . 

The analysis of the metallic line spectrum requires the best 
possible knowledge of atomic line parameters, log gf values in par- 
ticular. In this work we only use lines of Ca, Ti, Cr, and Fe for which 
experimental atomic parameters are available (except for Cm, as 
clarified later in the text). Atomic parameters were extracted from 
the VALD database and from other recent publications, see refer- 
ences in Table|4] 

Data for neutral an d ionised Ca lines were validated with 
non-LTE calculations by Masho nkina et ail d2007l) . For Til, Tin, 
Cri and Fei lines, accurate laboratory data (lifetimes and branch- 
ing ratio) are available. For lines of ionised iron the oscillator 
strengths, selected fr om VALD, were p r oduced from laboratory 
data, as explained in iRvabchikova et alj £l999). All correspond- 
ing references are given in Table [4] Laboratory data for optical 
Cm lines are scarce, therefore we took into account two differ- 
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ent sets of calculated data, t he one from semi-empirica l orthog- 
onal operator calculations by iRaassen & Uvlingj ( fl993) . and the 
one from the latest calculations by R. Kurucz. The lifetimes cal- 
culated by the two groups agree very well and are very c l ose to 
rece nt laboratory lifetim e measurements by iNilsson et alj ( 2006) 
and iGurell et al] J2010h . A small difference between theoretical 
and experimental lifetimes, converted to a difference in oscilla- 
tor strength, corresponds to a log gf- value uncertainty no larger 
than ±0.03 dex. However, the two sets of calculated oscillator 
strengths, for the lines used in our analysis, differ by 0.2 dex, as 
a result of different branching factors, as the lifetimes are prac- 
tically identical. Fr o m two theoretical sets of oscillator strengths, 
IRaassen & Uvlingj |l998)'s data provide us with a smaller stan- 
dard deviation from the mean Cm abundance, while Kurucz's data 
give a mean abundance closer to that derived from the Cri lines. 
Clearly, an extensive laboratory analysis of Crn lines in the opti- 
cal region is needed for the interpretation of the Cr abundance in 
atmospheres of cool and hot stars. 

LTE line abundances w ere based on eq uivalent widths, anal- 
ysed with a mo dified version (Tsymb all 19961) of the WIDTH9 code 
kuruczll 19931) . For blended lines and lines situated in the wings of 
the hydrogen lines we derived the line abundance performing syn- 
thetic spectrum calculations with the SYNTH3 code, tuning y macro 
line by line (see Sect. l4.2.2l . A line-by-line abundance list with the 
equivalent width measurements, adopted oscillator strengths, and 
their sources is given in Table [4] (see the online material for the 
complete version of the table). Table [4] also gives equivalent width 
measurements for the lines which we measured via spectral synthe- 
sis. In these cases, equivalent widths were tuned to match the abun- 
dance obtained with spectral synthesis. Our analysis shows that it 
is practically impossible to get a unique value of the microturbulent 
velocity for all considered species, therefore we derive a value that 
satisfies all the data and still provides a small scatter. 

Figures [4] and [5] show the correlations of Fei and Fen abun- 
dance with equivalent width (left panels) and with excitation po- 
tential (right panels), respectively for HD 321 15 and HD 37594. In 
each Figure, we used our final adopted fundamental parameters for 
the top panels and B08's fundamental parameters for the bottom 
panels. B08 derived the effective temperature by imposing the ex- 
citation equilibrium for the Fei lines only, the surface gravity by 
imposing the ionisation equilibrium for Fe i only, and the microtur- 
bulence velocity by removing the correlation between abundance 
and equivalent widths for Fe i lines only. 

With our model parameters for HD 321 15, we get a small pos- 
itive correlation for Fei abundance, and a small anti-correlation 
for Fe ii abundance, with the excitation potential. Our adopted T c g 
thus accommodates both Fei and Fen. With our model parame- 
ters for HD 37594, we get a small positive correlation for both Fe i 
and Fe n abundance with the excitation potential. These correlations 
would change slightly by adopting a different set of Fe lines, even 
by adding or removing a few lines, making the parameter deter- 
mination based on these correlations rather sensitive to systematic 
effects introduced by the line selection. One example is the high ex- 
citation energy Fen line, shown in Figs.|4]and[5] which has a much 
larger influence on the Fe n excitation equilibrium, compared to the 
other Fe n lines. 

The model parameters obtained by B08 lead to almost perfect 
equilibria for Fe i, as requested by B08 analysis method; in the case 
of HD 321 15, B08's r eff is anyway so high that even the excitation 
equilibrium for Fe i is not reached, although it is imposed by their 
analysis method. 

The equilibria for Fei are reached better with B08's parame- 
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Figure 4. Individual abundances for Fe i (black crosses) and Fe n (red as- 
terisks) lines versus the equivalent width (left panels) and the excitation 
energy of the lower level (right panels) for HD 32115, with different model 
atmospheres. The top panels are for our adopted fundamental parameters, 
while the bottom panels are for B08's fundamental parameters. The lin- 
ear fit to the Fe i (black solid line) and Fe n (red dashed line) data are also 
shown. The black cross and red asterisk at negative equivalent widths show 
respectively the mean Fe i and Fe n abundance with their standard deviation. 
The black and red arrows show the direction of the non-LTE corrections for 
Fei (upwards/higher abundance) and Fen (downwards/lower abundance), 
respectively. 
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Figure 5. Same as in Fig.g] but for HD 37594. 



ters, compared to ours, as it is the only T c ff and f m j c indicator they 
used, but it leads to a set of parameters which does not fit all the 
other parameter indicators, first of all Fe n. 

The average abundances, shown in Figs. [4] and [5] show that 
we obtain the ionisation equilibrium for Fe, within the error bars. 
We notice also that our average Fen abundance is systematically 
higher by a few 0.01 dex, compared to Fe i, in agreement with non- 
LTE calculations by Mashonkinal d201ll) (the direction of the Fe 
non-LTE corrections is shown in Figs. [4] and [5). Adopting B08's 
parameters, we obtain a systematically lower Fen abundance, com- 
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Table 2. Lines used for the parameter determination. Wavelengths and excitation potentials are taken from the VALD database. The adopted log gf values 
are taken from different sources which are listed in the last column. "S" in the equivalent with column denotes line abundances det ermined by fitt i ng the 
observed line profile, with the equivalent width determined from the line abundance. The log gf values by Blackwell et al. 1 1980) and Baschek et alj jl970h 
were corrected by +0.2 and +0.16, respectively. For Cm the results obtained with two different sets of log gf values are given (see the Sect. 13.23V For 
each ion, the last line gives the average abundance and the standard deviation, with the number of lines in parenthesis. The entire table can be viewed in the 
electronic version of the Journal. 
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pared to Fe i, and the non-LTE corrections by Masho nkiru] to 111) 
would then worsen the ionisation equilibrium, indicating that B08's 
sets of parameters is inappropriate for these stars. 

Figure [6] shows the line abundance versus the equivalent 
widths for Ca, Ti, and Cr, in HD32115 and HD 37594, adopting 
our final fundamental parameters. As for Fe, these elements display 
a variety of small correlations and on average our fundamental pa- 
rameters are the ones which suit at best all considered ions. The 
ionisation equilibrium is reached within the error bars for all ele- 
ments considered here, except Cr, for which the equilibrium is ob- 
tained adopting Kurucz's oscillator strengths, as shown in Table|4] 
Also with laboratory data for Cr lines, we expect an improvement 
in the Cr ionisation equilibrium. 



4 DISCUSSION 

For HD32115 our final adopted set of parameters is: r eff = 
7250±100K, logg = 4.2±0.1 and u mlc = 2.5±0.2kms~\ For 
HD 37594 we derived T cff = 7150±100K, logg = 4.2±0.1 and 
L> mic = 2.6±0.2kms~'. The measured projected rotational veloc- 
ity (usini) is 8.3 and 17.0 km s -1 , respectively for HD32115 and 
HD 37594. The observed line profiles required also a substan- 
tial macroturbulence velocity (v macTO ), generally between 8 and 
lOkms 1 for both stars (see Sect. |4.2l2l for more details). These 
high values of i/ macro are i n line with an extrapola tion of the iWro- 
r eff co rrelations given bv lValenti & Fisherl ([2005 ) and lBruntt et al.l 

Our fundamental parameters are not perfect, by definition, but 
they provide, within the error bars, the best description of all avail- 
able observables: photometric colors, hydrogen and metallic line 
profiles. 
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Figure 6. Individual abundances for Ca (top panels), Ti (middle panels), 
and Cr (bottom panels) lines versus the equivalent width measured for 
HD32115 (left panels) and HD 37594 (right panels). In each panel, the 
black crosses indicate the lines of neutral elements, while the red asterisks 
indicate the lines of singly ionised elements. All abundances are derived as- 
suming our adopted fundamental parameters. The linear fits to the data are 
also shown. The black cross and red asterisk at negative equivalent widths 
show the average abundance respectively for the neutral and singly ionised 
elements, with their standard deviation. The poor Cr ionisation equilibrium 
is most likely due to the poor quality line data. 



4.1 Model fluxes and observed photometry 

Other methods which should always be used to check the funda- 
mental parameters obtained by spectroscopic means are: i) com- 
parison of synthetic and observed optical colors; ii) comparison be- 
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tween photometry, calibrated in flux, and synthetic spectral energy 
distribution (SED). We performed these comparisons for HD 321 15 
and HD 37594, aiming also to estimate the precision which can be 
reached with these two methods. 



4.1.1 Synthetic colors 

Table [3] summarises the comparison between observed and theo- 
retical color-indexes for different photometric systems for the two 
stars, adopting the fundamental parameters obtained in this work 
and by B08. We highlight the better agreeing theoretical index in 
each case. Overall out values provide better mat ches. Al l color s 
were calculated using modified computer codes by |Kuruc3 d 19931) . 
which take into account transmission curves of individual photo- 
metric filters, mirror reflectivity and a photomultipli er response 
funct ion. In contrast to Kurucz's procedures (Rel vea & Kurucj 
119781) which are based on the low resolution theoretical fluxes, our 
synthetic colors are computed from energy distributions sampled 
with a fine wavelength step, so integration errors are expected to be 
small. 

To fully understand Table [3] it is first necessary to estab- 
lish more in general how well synthetic colors can reproduce 
the observations. This can be done by comparing synthetic and 
observed colors for well known stars. To check the quality of 
synthetic and observed colors over a large temperature range, 
we examine three "st andard" stars with d ifferent temperatures: 
Procyon (r. ff =6530K;lFuhrmann et al.ll 19971) , Vega ( r eff =9550 K; 
Castelli & Kuruczll 19941) . and 21 Peg (T eff = 10400 K; iFossati et al.l 
20091) . 



Figure|7]shows a comparison between synthetic and observed 
colors for the three reference stars, plus HD 321 15 and HD 37594. 
With a few exceptions, there is general good agreement between 
synthetic and observed colors throughout the temperature region 
explored here. For 21 Peg, there is difference of 0.04 mag between 
synthetic and observed U — B Johnson color, V\ -B and G-B Geneva 
colors. This difference is most likely due to incorrect photometry 
in one (or more) photometric band, as there is a perfect agreement 
between synthetic spectral energy distribution and spectrophotom- 
etry in the whole waveleng th region between t he near-UV and the 
near infrared (see Fig. 4 by Fossati et al. 2009). On the other hand 
we do not have an explanation for the difference of 0.03 mag ob- 
tained between synthetic and observed b — y Stromgren color for 
Procyon, although an error in the photometry is always possible, 
even for such a bright star. Besides these two exceptions, there is 
general agreement between synthetic and observed colors at a 0.01- 
0.02 mag level, which represents then the typical precision that one 
can expect for the comparison between observed and synthetic col- 
ors. 

Figure |7] shows that for HD 32115 and HD 37594 the differ- 
ence between synthetic and observed colors is within the typical 
precision obtained for the "reference" stars, confirming the quality 
of our fundamental parameters. Figure [8] shows that a comparison 
between observed and synthetic colors calculated with B08's stellar 
parameters clearly demonstrate the inadequacy of the B08 funda- 
mental parameters. 



4.1.2 Spectral energy distribution 

For a complete self-consistent analysis of any star, one should re- 
produce the observed spectral energy distribution with the adopted 
parameters for a model atmosphere. For HD 321 15 and HD 37594 
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Figure 7. Comparison between observed and synthetic colors for a set of 
reference stars, plus HD 321 15 (bright blue triangles) and HD 37594 (violet 
squares). The reference stars are: Vega (black crosses), Procyon (red aster- 
isks), and 21 Peg (blue diamonds). Photometric colors are labelled on the 
X-axis. For all stars we assumed zero interstellar reddening. 
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Figure 8. Comparison between observed and synthetic colors for HD 321 15 
(shifted downwards of 0.1) and HD 37594 (shifted upwards of 0.1). The 
colors calculated with our final set of parameters and B()8's parameters are 
shown by black crosses and red asterisks, respectively. Photometric colors 
are labelled on the xaxis. For both stars we assumed zero interstellar red- 
dening. 



no spectrophot ometry is available, th erefore w e converted the a vail- 
able Johnson ( Johnson et al ] |l966l) . Geneva llRufeneri[T988l) and 
2MASS dZacharias et alj2005l) photometry into physical units and 
compared them with the model fluxes calculated with the final 
set of parameters derived for the two stars. In the case of the 
Geneva colors, we assumed that Geneva V G index is close to John- 
son Vj, making then possible to recover the Geneva B index and 
thus U, Bu Bi, Vi, and G, and to t ransform them to absolut e 
fluxes using the calibrations given by iRufener & Nicoletl dl988l) . 
For the 2MASS (JHK) photometry we employed calibrations by 
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Table 3. Observed and calculated photometric parameters of HD 321 15 and HD 37594. The values in brackets give the error bars of observations. The better 
agreeing theoretical index in each case is highlighted. 
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van der Bliek et al while for the Joh nson photometry we 

employed calibrations bv lBessell et al. (1998). 

Figure [9] shows the comparison between the observed pho- 
tometry (in physical units) and the model fluxes calculated with 
the adopted atmospheric parameters for HD32115 and HD 37594. 
There is a good agreement between observed and synthetic fluxes 
in the visible and infrared regions, providing a further confirmation 
of our adopted stellar parameters for the two stars. 

Precise H IPPARCOS parallaxes are available for both stars 
dvan Le euwen 2007). This allowed us also to estimate their radii: 
1.52±0.04J? o for HD32115 and 1.36±0.02# o for HD37594. 
The radius we derived for HD32115 i s in very good agree- 
ment with what previous l y obt ained by iFekel et alT l feOOrjl) and 
lAUende Prieto & Lambert. Jl999h . who measured the stellar radius 
by means of stellar evolution calculations. This comparison in- 
creases the confidence on our results. 
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4.2 Classic vs. semi-automatic procedures 

The most evident difference between the analysis performed by 
B08 and our analysis (including HD 49933 by R09) is that B08 
did not add any L> macro to the line broadening. Since B08 used 
line profile fitting, instead of equivalent widths, the lack of the 
additional broadening renders their whole analysis questionable. 
Bik maev et"al] j2002h also did not use any v mllao in their analysis, 
but their spectra were of mid-low resolution and therefore vsini 
was enough to fit the lines and their abundance analysis was based 
mainly on equivalent widths, independent of u mxT0 . At the resolu- 
tion of CORALIE and TS (R=50 000 and 60 000, respectively), line 
profiles clearly require an additional f macro broadening. 

Figure [TOl compares the observed HD32115 line profiles of 
two Fen lines, one weak and one strong, with three synthetic spec- 
tra calculated with our adopted parameters and three different sets 



Figure 10. Observed HD32115 line profiles (black solid line) of two Fen 
lines at /15132A (left panel) and /15197A (right panel) in comparison 
with three synthetic spectra calculated with three different sets of broad- 
ening parameters: Msinr'=8.7kms~' and f macro =8.0kms -1 (red dashed 
line); vsini=8.7kms 1 and (J macro =9. 8kms-' ( green dotted-dashed line); 
v sin f'=9.0 km s and v macro 

=0 km s 1 (blue dotted line). The first set of pa- 
rameters fits the weak lines best, the second the strong lines, while the third 
is the one adopted by BOS. For all three synthetic spectra we used our final 
set of fundamental parameters and abundances. 



of broadening parameters. The three broadening parameter sets are 
suited to weak lines, strong lines, and one as given by B08, respec- 
tively. 

The difference in u maclo between weak and strong lines is 
found systematically for several lines of the same ion, and for 
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Figure 9. Comparison between LLmodels theoretical fluxes (full red lines) calculated with the fundamental parameters derived for HD32115 (upper panel) 
and HD 37594 (lower panel), with Johnson, Geneva (respectively open diamonds and squares in the left panels) and 2MASS (open triangles in the right panels) 
photometry converted to physical units. 



different ions, excluding the possibility that this is due to an er- 
ror in v mic , v sin i, or in the damping constants (homogeneous cal- 
culati ons bv lBarklem et alJboOCt fBarklem & Aspelund- Johansson! 
2005). This systematic diffe rence in i/ maC ro betwe en weak and 
strong lines was found also bv lFuhrmann et all d 19971) in their anal- 
ysis of Procyon. 

The profiles calculated using the broadening by B08 show too 
deep cores and too narrow wings. This implies that B08 systemati- 
cally obtained erroneously high r c ff for the three stars (HD 49933, 
HD 321 15, and HD 37594) to compensate for the deeper line cores, 
as a result of neglecting y macl0 and fitting line c ores to derive the 
abundances. This is confirmed by the fact that iBrurit tl j2009h re- 
analysed HD 49933, adopting a i> maclo of 2kms~ 1 , and obtained a 
much lower temperature than given in B08. It is still not clear to us 
where this particular value of lw io came from, as R09 measured 
a f m a C ro of 5.2±0.5 km s -1 , in agreement with the calibration by 
IValenti & Fisiien (2005). Figure [TOl demonstrates that even adopt- 
ing all the parameters and abundances given by B08, it is impossi- 



ble to simultaneously fit the line profiles of weak and strong lines, 
as the wings will always be too narrow and the higher T e n cannot 
compensate for it. 

One of our main goals was to find and analyse possible 
sources of discrepancies between the results obtained with a clas- 
sical method of analysis (e.g. this work and R09) and with a semi- 

I — — I 1 1 1 71. 

automatic procedure (e.g. Bruntt et al. 2008; Bruntt 2009). We have 

identified two possible sources of discrepancy. The first arises from 

semi-automatic procedures often taking into account only Fei lines 

to derive r c(f and u mic , while the second arises from the nature of 

the line profiles for these stars. 



4.2. 1 Is Fe I alone good enough for a precise parameter 
determination ? 

The upper-left panel of Fig. [4] shows the Fe i and Fe n line abun- 
dance as a function of equivalent width, used for the determination 
of the microturbulence velocity. It is clear that a perfect equilibrium 
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between the line abundance and equivalent width is not reached, 
and not simultaneously reachable, as the two ions show opposite 
correlations: -5.524 ± 2.407 and 6.402 ± 3.676, in units of 10~ 4 , re- 
spectively for Fei and Fen. Something similar is found for the ions 
of other species, as both Ca and Cr show a positive correlation, 
while Ti shows a negative correlation (see Fig.[6j. 

The Vaac value we adopted considers the effects on all the ions 
we took into account, whereas the v mic value derived from the Fe i 
lines alone, worsens the situation for Fen and Ti, re sulting in a 
systematic underestimation of f mi ,.. IBruntt et al.l {2010) presented a 
calibration of L> m j C as a function of T c g for late-type stars, based on 
the analysis of Fe i lines alone, therefore this calibration could also 
underestimate v mic . 

Using Fe i lines alone for the determination of v mic means that 
the average Fe n abundance depends strongly on which lines have 
been selected. If the selected Fen lines are predominantly weak, 
the Fe n abundance will be artificially low. This introduces a bias 
in the determination of the surface gravity, as most semi-automatic 
procedures employ the Fe ionisation equilibrium alone to measure 
logg. The use of predominant medium-to-weak Fen lines, with the 
adoption of Fe i lines alone to determine L» m j c consequently leads to 
a sistematically erroneously low Fe n abundance, which then causes 
an erroneously high log g to be inferred. This of course, then affects 
the value of r e ff- 

Another problem connected to the use of Fe i lines as the sole 
temperature indicator is shown in the top-right panel of Fig. [4] 
The effective temperature which best fits all available observables 
does not completely remove the correlation between abundance 
and excitation potential. For HD32115, by adopting only Fei as 
a r e ff and u m ; c indicator, and only Fe for the ionisation equi- 
librium, we derived r cff =7400K, logg=4.2, and u mic =2.4kms~ 1 . 
Similarly, for HD 37594 we obtained r cff =7300K, logg=4.1, and 
t>mic=2.5 km s~ . These temperatures are too high and cannot fit the 
hydrogen lines or the photometry, with photometric colors reveal- 
ing the discrepancies. 

The adoption of Fe i alone as T eff and v mic indicator might also 
lead to the presence of systematic errors in the fundamental pa- 
rameters arising from non-LTE effects. This is particularly true for 
cool stars, where non-LTE effects are larger for Fe i than Fe n, with 
obvious consequences for i.e. logg. Non-LTE effects are generally 
stronger for stronger lines, implying that systematic errors might 
be introduced in the determination of u m ; c , consequently affecting 
r c ff as well. 

All problems described here can be solved by including other 
ions in the process of parameter determination. Although fewer 
lines of ions other than Fei are usually available, their inclusion 
in the parameter determination procedure would significantly alle- 
viate the systematic errors introduced by the use of only Fei lines. 

lAsplund et all d2005h shows differences between 3D and ID 
LTE abundances for several ions as a function of the excitation po- 
tential for the Sun. For Fei, the ID models, compared to the 3D 
models, lead to a deviation as large as 0.2 dex for the low excitation 
lines, which decreases with increasing excitation potential. Clearly, 
the use of ID models introduces a strong bias towards high temper- 
atures. This effect is present, with slightly different strengths, for 
many other ions. It would be extremely valuable to study this effect 
for stars hotter than the Sun, where the registered high v mmm values 
suggest that hydrodynamical effects might be even stronger than in 
the Sun. 



4.2.2 The effect of macroturbulence 

IValenti & Fished d2005l) and iBruntt et al.l d2010h showed that for 
late-type stars the macroturbulence velocity (f macro ) increases with 
increasing T eS . For HD 321 15 and HD 37594, i> macro reaches rather 
high values, close to lOkms -1 . During the analysis of the spectra 
of the two stars, we noticed a systematic difference between the 
line profiles of strong and weak lines, as strong lines require a sys- 
tematically higher iWm compared to weak lines, by ~2km s~ , as 
shown in Fig.[l0] We believe that this difference is due to unmod- 
elled depth-dependent velocities in the atmosphere which become 
evident for stars with relatively high u macro and low v sin i. 

This is not a concern if the abundances are derived from equiv- 
alent widths, as they are independent of macroturbulence velocity. 
On the other hand, if the abundances are obtained from line profile 
fitting, where the abundance is the only free parameter, a system- 
atic abundance difference is introduced between strong and weak 
lines, affecting the determination of the fundamental parameters. 

Semi-automatic abundance analysis procedures first set the 
broadening parameters to the entire available spectrum and then 
measure the line abundance from line profile fitting. If strong and 
weak lines require different broadening, opposing biases are in- 
troduced in the abundances derived from weak and strong lines. 
We estimated that the abundances obtained by fitting the profile of 
strong lines, adopting a v mact0 typical for weak lines, are systemati- 
cally lower by 0. 10-0. 15 dex, compared to the abundances obtained 
with a u raacro appropriate to strong lines. Consequently the abun- 
dances obtained from weak and strong lines will be systematically 
over- and under-estimated, respectively. 

The systematic effect described above might be temperature 
and u sin i dependent as the difference in f macr0 between weak and 
strong lines could increase with increasing r cff , and be more evi- 
dent with decreasing v sin i. In fact, the worst situation may be for 
slowly rotating "cool" early-type stars, such as HD 321 15. 

The problem described here can be solved using equivalent 
widths, instead of line profile fitting, as the equivalent widths are 
broadening independent. Alternatively u macro could be treated as 
a further free parameter in the line profile fitting, although thor- 
ough tests on the uniqueness of the derived line abundance and 
Vmmo should be performed. The adoption of 3D atmosphere mod- 
els would also partly eliminate this problem, but this is not a viable 
solution for an abundance analysis of a large sample of stars. 



5 CONCLUSION 

On the basis of high resolution, high signal-to-noise ratio spec- 
tra, taken at the McDonald observatory, we carried out a pre- 
cise parameter determi nation of two late A -type stars, HD32115 
and HD 37594, which IBruntt et all |2008) adopted as reference 
stars for their abundance analysis of yDor stars. IBruntt et al.l 
(2008) analysed these stars with a semi-automatic procedure and 
the ir results stron g ly dis agreed with those previously obtained 
by iBikmaev et alj | |2002|) . by means of a "classical" analysis. 
This discrepancy, togeth er with that previously highlighted by 
iRvabchikova et all d2009f) for HD 49933, prompted us to reanalyse 
HD 321 15 and HD 37594 to look for the origin of the discrepancies. 

We derived the fundamental parameters making use of all the 
available observables: multicolor photometry, pressure-sensitive 
magnesium lines, metallic lines and profiles of hydrogen Balmer 
lines. For HD32115 our final adopted set of parameters is: T tf! 
= 7250±100K, logg = 4.2+0.1 and v mic = 2.5+0.2 km s" 1 . For 
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HD 37594 we adopted: T cff = 7150±100K, logg = 4.2+0.1 and 
v mic = 2.6±0.2km s~ . We confirmed our final set of parameters by 
comparing flux calibrated photometry with synthetic spectral en- 
ergy distributions, and by comparing observed and synthetic pho- 
tometric colors. Our prefered fundamental parameters fit, within 
the error bars, all availa ble observables . They are also in agree- 
ment with th e results by Bikmaev et alj ( 120021) . but disagree with 
the results bv lBruntt et alj J2008). For HD 321 1 5, our set of param- 
eter s agrees with that previou s ly obt ained by iFekel et alj J2006h 
and lAllende Prieto & Lambert! {1999). In addition, the radius of 
HD 321 15, we derived from the analysis of the spectral energy dis- 
tributi on, is in perfect agreement with that d erived by IFekel et alj 
(2006) and lAllende Prieto & Lambertl dl999r) by means of stellar 
evolution calculations. 

We co mpared our McDonald spectra with those used by 
iBruntt et alj ( 120081) and lBikmaev et all fe002h , concluding that the 
differences between the three set of spectra is well within the given 
S/N. The quality of the data is not the origin of the discrepancies. 

To fit the line profiles of the two stars, we had to adopt rather 
large u maclo values, between 8 a nd lOkms -1 . Although the spectra 
analysed by IBruntt et alj d2008l) had enough resolution to require 
the need of u macro to fit line profiles, they adopted only rotational 
broadening (v sin z) in their spectral synthesis. As a consequence, 
by deriving the line abundance by fitting synthetic line profiles to 
the observed ones, discarding any i/ macro broadening, they intro- 
duced a bias in their analysis, which we believe led them to an 
erroneous set of fundamental parameters. 

We have demonstrated that the determination of T eff and log g 
using only the Fei excitation equilibrium and the Fe ionisation 
equilibrium leads to a systematic higher compared to that sug- 
gested by all other indicators. We also believe that this effect might 
be temperature dependent. 

Several automatic and semi-automatic procedures use the 
Fei excitation equilibrium and the Fe ionisation equilibrium as 
only/primary indicators for stellar parameter determination. In this 
work we show that these procedures do not always provide the cor- 
rect set of fundamental parameters and their results need to be cau- 
tiously taken when performing further analysis, such as modelling 
of the asteroseismic frequencies or the characterisation of transiting 
exoplanets. 
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Table 4: Lines used for the parameter determination. Wavelengths and excita- 
tion potentials are taken from the VALD database. The adopted logg/ values are 
taken from different sources which are listed in the last column. "S" in the equiv- 
alent with column denotes line abundances determined by fitting the observed 
line profile, wit h the equivalent width dete rmined from the line abundance. The 
log gf values bv lBlackwell et al! | imdt) and lBaschek etail dl970h were corrected 
by +0.2 and +0.16, respectively. For Cm the results obtained with two different 
sets of logg/ values are given (see the Sect. |3.23t . For each ion, the last line 
gives the average abundance and the standard deviation, with the number of lines 
in parenthesis. 



Element 






HD 37594 


HD32115 




Wavelength 


^excit 




EQW 


abundance 


EQW 


abundance 


Ref log gf 


A 


eV 




mA 


dex 


mA 


dex 




Cai 
















4425.4370 


1.8790 


-0.358 


S 100.0 


-6.00 


112.71 


-5.70 


SN 


4435.6790 


1.8860 


-0.517 


S90.0 


-6.00 






SN 


4526.9280 


2.7090 


-0.548 






54.80 


-5.78 


SR 


4578.5510 


2.5210 


-0.697 


43.00 


-5.99 


54.75 


-5.75 


SR+Sm 


4685.2680 


2.9330 


-0.879 


21.46 


-5.90 


27.60 


-5.70 


S 


5261.7040 


2.5210 


-0.579 


S46.0 


-6.08 






SR+Sm 


5512.9800 


2.9330 


-0.464 


39.77 


-6.00 


47.26 


-5.83 


SR+Sm 


5581.9650 


2.5230 


-0.555 


55.26 


-5.97 


64.10 


-5.78 


SR+Sm 


5588.7490 


2.5260 


0.358 


125.68 


-5.84 






SR+Sm 


5590.1140 


2.5210 


-0.571 


57.83 


-5.92 


64.85 


-5.75 


SR+Sm 


5601.2770 


2.5260 


-0.523 


56.75 


-5.98 


73.14 


-5.68 


SR+Sm 


5857.4510 


2.9330 


0.240 


101.96 


-5.81 






S+Sm 


5867.5620 


2.9330 


-1.570 


4.87 


-5.94 


7.73 


-5.66 


S 


6122.2170 


1.8860 


-0.316 


123.55 


-5.80 






SN 


6161.2970 


2.5230 


-1.266 


13.42 


-6.09 


23.49 


-5.74 


SR+Sm 


6162.1730 


1.8990 


-0.090 


144.67 


-5.74 






SN 


6166.4390 


2.5210 


-1.142 


24.02 


-5.92 


30.20 


-5.73 


SR+Sm 


6169.0420 


2.5230 


-0.797 


32.57 


-6.10 


52.80 


-5.71 


SR+Sm 


6169.5630 


2.5260 


-0.478 


72.66 


-5.83 


74.77 


-5.73 


SR+Sm 


6471.6620 


2.5260 


-0.686 


45.97 


-5.99 


64.32 


-5.66 


SR+Sm 


6493.7810 


2.5210 


-0.109 


99.99 


-5.80 






SR+Sm 


6499.6500 


2.5230 


-0.818 


41.29 


-5.93 






SR+Sm 


6717.6810 


2.7090 


-0.524 


56.24 


-5.87 






SR 


7148.1500 


2.7090 


0.137 


107.26 


-5.80 


117.58 


-5.58 


SR 


7202.2000 


2.7090 


-0.262 


76.54 


-5.85 






SR 


7326.1450 


2.9330 


-0.208 


67.13 


-5.87 


85.61 


-5.65 


S 




Average 




-5.92+0.10(25) 


-5.71+0.06 (16) 




Can 
















5001.4790 


7.5050 


-0.507 


S 36.00 


-5.82 


48.88 


-5.63 


TB 


5019.9710 


7.5150 


-0.247 


S 45.50 


-5.92 






TB 


5021.1380 


7.5150 


-1.207 


S7.60 


-5.92 


8.00 


-5.91 


TB 


5285.2660 


7.5050 


-1.147 


S 11.10 


-5.80 


S 16.40 


-5.62 


TB 


5339.1880 


8.4380 


-0.079 


S 25.50 


-5.79 


S 40.00 


-5.60 


TB 


6456.8750 


8.4380 


0.412 


S 61.00 


-5.64 


80.00 


-5.58 


TB 


8201.7220 


7.5050 


0.368 


S 100.00 


-5.72 


125.03 


-5.50 


TB 


8248.7960 


7.5150 


0.556 


S 163.00 


-5.57 






TB 


8254.7210 


7.5150 


-0.398 


S 57.00 


-5.72 


68.79 


-5.60 


TB 




Average 




-5.77+0.12(9) 


-5.64+0.13(7) 




Til 

4453.6990 


1.8730 


-0.010 






9.30 


-7.07 


MFW 


4548.7630 


0.8260 


-0.354 


17.07 


-7.34 


20.80 


-7.14 


MFW 


4617.2690 


1.7490 


0.389 


13.00 


-7.50 


19.90 


-7.20 


MFW 


4758.1180 


2.2490 


0.425 






10.80 


-7.16 


MFW 


4759.2700 


2.2560 


0.514 






13.38 


-7.14 


MFW 


4981.7310 


0.8480 


0.504 


60.55 


-7.44 


68.20 


-7.24 


MFW 
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Element 






HD 37594 


HD32115 




Wavelength 




logg/ 


EQW 


abundance 


EQW 


abundance 


Ref log gf 


A 


eV 




mA 


dex 


mA 


dex 




4999.5030 


0.8260 


0.250 


49.47 


-7.36 


51.70 


-7.23 


MFW 


5192.9690 


0.0210 


-1.006 


22.40 


-7.23 


21.30 


-7.15 


MFW 


5210.3850 


0.0480 


-0.884 


20.84 


-7.36 


23.10 


-7.21 


MFW 


6261.0980 


1.4300 


-0.479 


7.67 


-7.19 


8.20 


-7.07 


MFW 




Average 




-7.35+0.11 (7) 


-7.16+0.06(10) 




Tin 
4394.0510 


1.2210 


-1.770 






93.50 


-7.09 


BHN 


4395.8389 


1.2430 


-1.970 






74.00 


-7.15 


BHN 


4411.9250 


1.2240 


-2.520 


33.86 


-7.26 


39.00 


-7.11 


PTP 


4417.7136 


1.1650 


-1.190 


133.34 


-7.14 


138.58 


-6.95 


PTP 


4418.3300 


1.2370 


-1.990 


70.91 


-7.25 


75.50 


-7.12 


BHN 


4421.9380 


2.0610 


-1.660 


50.88 


-7.19 


59.00 


-7.02 


PTP 


4444.5545 


1.1160 


-2.210 


63.98 


-7.22 


66.00 


-7.12 


BHN 


4470.8532 


1.1650 


-2.020 






76.16 


-7.14 


PTP 


4568.3140 


1.2240 


-2.940 






21.63 


-7.03 


PTP 


4609.2640 


1.1800 


-3.430 


8.31 


-7.10 


7.07 


-7.11 


BHN 


4636.3200 


1.1650 


-3.230 


7.18 


-7.38 


12.35 


-7.07 


BHN 


4657.2004 


1.2430 


-2.320 






51.03 


-7.13 


BHN 


4708.6621 


1.2370 


-2.370 


40.00 


-7.32 


48.70 


-7.12 


BHN 


4779.9850 


2.0480 


-1.260 






81.50 


-7.15 


RHL 


5005.1570 


1.5660 


-2.730 


13.78 


-7.28 


17.44 


-7.10 


BHN 


5185.9018 


1.8930 


-1.490 


S69.0 


-7.30 


82.00 


-7.07 


PTP 


5252.0190 


2.5900 


-1.960 


13.79 


-7.24 


14.86 


-7.16 


PTP 


5336.7710 


1.5820 


-1.630 






87.09 


-7.12 


BHN 


5381.0150 


1.5660 


-1.970 


56.00 


-7.26 


64.86 


-7.09 


BHN 


7214.7160 


2.5900 


-1.740 


19.60 


-7.35 


26.44 


-7.01 


MFW 




Average 




-7.25±0.08 (13) 


-7.09+0.05 (20) 




Cri 
















4274.7970 


0.0000 


-0.231 


S 142.0 


-6.72 


S 153.0 


-6.37 


MFW 


4545.9530 


0.9410 


-1.370 


22.29 


-6.84 


S 34.0 


-6.51 


SLS 


4492.3050 


3.3750 


-0.392 






9.32 


-6.29 


MFW 


4616.1240 


0.9830 


-1.190 


28.18 


-6.87 


41.63 


-6.53 


SLS 


4646.1620 


1.0300 


-0.740 






68.94 


-6.56 


SLS 


4651.2840 


0.9830 


-1.460 


18.32 


-6.83 


27.65 


-6.51 


SLS 


4652.1570 


1.0040 


-1.040 


41.90 


-6.76 


53.42 


-6.50 


SLS 


4689.3570 


3.1250 


-0.400 


11.00 


-6.47 


15.33 


-6.24 


SLS 


4708.0130 


3.1680 


0.070 


19.50 


-6.62 


20.79 


-6.52 


SLS 


4718.4200 


3.1950 


0.240 


29.00 


-6.57 


30.42 


-6.46 


SLS 


4752.0870 


4.1860 


0.440 






16.68 


-6.24 


MFW 


4756.1120 


3.1040 


0.090 






34.43 


-6.31 


MFW 


4789.3350 


2.5440 


-0.330 






28.78 


-6.42 


SLS 


4922.2650 


3.1040 


0.380 






49.08 


-6.37 


SLS 


4936.3360 


3.1130 


-0.250 


11.50 


-6.61 


16.30 


-6.37 


SLS 


5206.0370 


0.9410 


0.020 






140.00 


-6.31 


SLS 


5247.5650 


0.9610 


-1.590 






26.00 


-6.46 


SLS 


5296.6910 


0.9830 


-1.360 


23.50 


-6.82 


32.78 


-6.55 


SLS 


5297.3770 


2.9000 


0.167 


21.91 


-6.89 


40.75 


-6.48 


MFW 


5348.3150 


1.0040 


-1.210 


31.00 


-6.82 


38.04 


-6.59 


SLS 


5409.7840 


1.0300 


-0.670 






84.09 


-6.46 


SLS 


5783.0630 


3.3230 


-0.500 






4.51 


-6.59 


MFW 


5787.9180 


3.3220 


-0.083 






19.30 


-6.32 


MFW 


6925.2720 


3.4490 


-0.330 






5.82 


-6.56 


MFW 


6978.3970 


3.4640 


0.142 






18.87 


-6.47 


MFW 


6979.7950 


3.4640 


-0.410 






6.42 


-6.43 


MFW 




Average 




-6.74+0.14(12) 


-6.44+0.11 (26) 




Cm 
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Element 






HD 37594 


HD32115 




Wavelength 


^excit 


logg/ 


EQW 


abundance 


EQW 


abundance 


Ref logg/ 


A 


eV 




mA 


dex 


mA 


dex 




4145.7810 


5.3190 


-1.106 






S42.5 


-6.10 


RU 


4252.6320 


3.8580 


-2.054 


S 33.0 


-6.44 


S51.0 


-6.12 


RU 


4275.5670 


3.8580 


-1.736 


S54.0 


-6.43 


S69.0 


-6.18 


RU 


4554.9880 


4.0710 


-1.491 


S56.0 


-6.50 


S77.0 


-6.16 


RU 


4588.1990 


4.0710 


-0.845 


S 107.0 


-6.39 


S 124.0 


-6.05 


RU 


4592.0490 


4.0740 


-1.473 


S66.0 


-6.37 


S 87.5 


-6.02 


RU 


4616.6290 


4.0720 


-1.576 


S56.0 


-6.45 


S70.0 


-6.18 


RU 


4618.8030 


4.0740 


-1.084 


S88.0 


-6.45 


S 111.0 


-6.03 


RU 


4634.0700 


4.0720 


-1.236 


S 81.0 


-6.40 


S 97.5 


-6.11 


RU 


5237.3290 


4.0730 


-1.350 


S71.0 


-6.45 


S 86.0 


-6.20 


RU 




Average 




-6.43+0.04 (9) 


-6.12±0 


07 (10) 




4145.7810 


5.3190 


-1.23 






S42.5 


-5.98 


K10 


4252.6320 


3.8580 


-1.99 


S 33.0 


-6.50 


S51.0 


-6.18 


K10 


4275.5670 


3.8580 


-1.67 


S 54.0 


-6.50 


S69.0 


-6.25 


K10 


4554.9880 


4.0710 


-1.28 


S56.0 


-6.71 


S77.0 


-6.37 


K10 


4588.1990 


4.0710 


-0.63 


S 107.0 


-6.61 


S 124.0 


-6.27 


K10 


4592.0490 


4.0740 


-1.22 


S66.0 


-6.64 


S 87.5 


-6.27 


K10 


4616.6290 


4.0720 


-1.36 


S 56.0 


-6.67 


S70.0 


-6.40 


K10 


4618.8030 


4.0740 


-0.83 


S 88.0 


-6.70 


S 111.0 


-6.28 


K10 


4634.0700 


4.0720 


-1.02 


S 81.0 


-6.62 


S97.5 


-6.33 


K10 


5237.3290 


4.0730 


-1.14 


S71.0 


-6.66 


S 86.0 


-6.41 


K10 




Average 




-6.62+0.08 (9) 


-6.31+0 


10(10) 




Fei 
4168.9416 


3.4170 


-1.650 






27.60 


-4.58 


FMW 


4189.5550 


3.6940 


-1.330 






25.69 


-4.73 


FMW 


4213.6474 


2.8450 


-1.290 


55.70 


-5.01 


71.11 


-4.69 


FMW 


4233.6020 


2.4820 


-0.604 






139.12 


-4.57 


FMW 


4248.2240 


3.0710 


-1.286 


48.90 


-4.94 


61.29 


-4.67 


BWL 


4250.1180 


2.4690 


-0.405 


146.00 


-4.79 


156.31 


-4.58 


FMW 


4266.9640 


2.7270 


-1.812 


36.67 


-4.87 


49.96 


-4.57 


BWL 


4267.8260 


3.1110 


-1.174 


56.61 


-4.91 






BWL 


4433.2170 


3.6540 


-0.700 


58.90 


-4.95 


80.53 


-4.58 


FMW 


4484.2190 


3.6020 


-0.864 


64.40 


-4.75 






BWL 


4485.6750 


3.6860 


-1.020 


35.83 


-4.95 


50.21 


-4.65 


FMW 


4602.0000 


1.6080 


-3.154 


16.70 


-4.86 


21.98 


-4.61 


FMW 


4602.9410 


1.4850 


-2.209 


73.70 


-4.94 


87.52 


-4.63 


BWL 


4630.1200 


2.2790 


-2.587 


18.00 


-4.86 


24.27 


-4.61 


BWL 


4635.8460 


2.8450 


-2.358 






15.50 


-4.63 


BWL 


4643.4630 


3.6540 


-1.147 


32.90 


-4.90 


41.41 


-4.68 


BWL 


4683.5597 


2.8310 


-2.319 






15.29 


-4.69 


BWL 


4690.1360 


3.6860 


-1.645 


14.00 


-4.83 






BWL 


4733.5910 


1.4850 


-2.988 






35.18 


-4.62 


BWL 


4736.7720 


3.2110 


-0.752 


85.92 


-4.90 


94.90 


-4.64 


BWL 


4745.8000 


3.6540 


-1.270 


27.93 


-4.88 


39.58 


-4.59 


BWL 


4789.6508 


3.5460 


-0.958 






58.77 


-4.70 


BWL 


4962.5719 


4.1780 


-1.182 


17.00 


-4.90 


24.49 


-4.57 


BWL 


4966.0870 


3.3320 


-0.871 


70.10 


-4.90 


92.31 


-4.48 


BWL 


4994.1290 


0.9150 


-3.080 


48.70 


-4.94 


55.75 


-4.69 


FMW 


5014.9410 


3.9430 


-0.303 


76.40 


-4.93 


87.60 


-4.64 


BWL 


5044.2100 


2.8510 


-2.038 


18.70 


-4.96 


28.09 


-4.66 


BWL,BK 


5049.8190 


2.2790 


-1.355 


84.50 


-5.03 


102.0 


-4.66 


BWL 


5054.6420 


3.6400 


-1.921 


5.00 


-5.09 


10.0 


-4.69 


BWL 


5083.3380 


0.9580 


-2.958 


48.20 


-4.99 


62.86 


-4.68 


FMW 


5090.7670 


4.2560 


-0.400 


46.49 


-4.99 


63.00 


-4.69 


FMW 


5127.3580 


0.9150 


-3.307 






41.22 


-4.68 


FMW 


5151.9100 


1.0110 


-3.322 


30.70 


-4.95 


40.87 


-4.59 


FMW 
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Element 






HD 37594 


HD32115 




Wavelength 


Afexcit 


logg/ 


EQW 


abundance 


EQW 


abundance 


Ref log gf 


A 


eV 




mA 


dex 


mA 


dex 




5194.9410 


1.5570 


-2.090 


71.06 


-5.01 


89.29 


-4.70 


FMW 


5198.7110 


2.2230 


-2.135 


35.30 


-5.02 


50.82 


-4.68 


FMW 


5232.9390 


2.9400 


-0.058 


158.60 


-4.85 


149.16 


-4.76 


BWL 


5236.2020 


4.1860 


-1.497 


5.00 


-5.09 


8.38 


-4.79 


BWL 


5242.4910 


3.6340 


-0.967 


50.40 


-4.83 


55.67 


-4.68 


BWL 


5269.5370 


0.8590 


-1.321 


167.67 


-4.86 


178.64 


-4.56 


FMW 


5281.7900 


3.0380 


-0.834 


80.31 


-5.03 


96.56 


-4.70 


BWL 


5283.6210 


3.2410 


-0.432 






111.91 


-4.71 


BKK 


5288.5247 


3.6940 


-1.508 


16.30 


-4.91 


23.18 


-4.65 


BWL 


5324.1780 


3.2110 


-0.103 


125.00 


-4.96 






BKK 


5364.8580 


4.4450 


0.228 


83.60 


-4.96 


101.52 


-4.66 


BWL 


5367.4660 


4.4150 


0.443 


96.60 


-5.04 


112.50 


-4.77 


BWL,BK 


5373.6980 


4.4730 


-0.860 


22.30 


-4.79 


33.26 


-4.51 


FMW 


5379.5740 


3.6940 


-1.514 


18.51 


-4.84 


23.33 


-4.64 


BWL 


5386.3330 


4.1540 


-1.770 


7.16 


-4.77 


8.98 


-4.51 


FMW 


5398.2770 


4.4450 


-0.670 


33.60 


-4.78 


44.39 


-4.54 


FMW 


5415.1920 


4.3860 


0.642 


118.70 


-4.95 


137.90 


-4.65 


BWL 


5434.5230 


1.0110 


-2.122 


103.70 


-5.02 


114.97 


-4.72 


FMW 


5473.9000 


4.1540 


-0.760 


38.99 


-4.88 


48.57 


-4.60 


FMW 


5497.5160 


1.0110 


-2.849 


68.04 


-4.81 






FMW 


5522.4460 


4.2090 


-1.550 


9.10 


-4.75 


13.71 


-4.49 


FMW 


5576.0888 


3.4300 


-1.000 


61.30 


-4.82 


75.26 


-4.55 


FMW 


5618.6310 


4.2090 


-1.276 


11.40 


-4.92 


19.03 


-4.60 


BWL 


5633.9460 


4.9910 


-0.270 


28.50 


-4.85 


41.42 


-4.57 


FMW 


5638.2620 


4.2200 


-0.870 


34.02 


-4.75 


42.79 


-4.53 


FMW 


5662.5160 


4.1780 


-0.573 






63.83 


-4.56 


BWL 


5679.0229 


4.6520 


-0.920 


20.20 


-4.65 


29.83 


-4.38 


FMW 


5701.5440 


2.5590 


-2.216 






31.85 


-4.66 


FMW 


5753.1210 


4.2600 


-0.688 






46.33 


-4.63 


BWL 


5816.3730 


4.5480 


-0.601 


29.12 


-4.93 


39.12 


-4.62 


FMW 


5852.2170 


4.5480 


-1.330 


9.10 


-4.72 


14.00 


-4.45 


FMW 


5905.6710 


4.6520 


-0.730 


19.42 


-4.90 


27.80 


-4.61 


FMW 


5916.2470 


2.4530 


-2.994 






11.41 


-4.50 


FMW 


5934.6530 


3.9280 


-1.170 


24.08 


-4.90 


35.43 


-4.58 


FMW 


6065.4820 


2.6080 


-1.530 


60.60 


-4.96 


75.56 


-4.67 


FMW 


6127.9060 


4.1430 


-1.399 


8.70 


-4.99 


16.49 


-4.61 


BWL 


6136.6150 


2.4530 


-1.400 


86.11 


-4.98 


91.65 


-4.69 


FMW 


6136.9930 


2.1980 


-2.950 






16.57 


-4.57 


FMW 


6137.6910 


2.5880 


-1.403 


75.81 


-5.00 


89.02 


-4.62 


FMW 


6151.6170 


2.1760 


-3.299 






6.51 


-4.68 


FMW 


6165.3600 


4.1430 


-1.474 


10.50 


-4.82 


14.46 


-4.60 


BWL 


6180.2030 


2.7270 


-2.586 


8.40 


-4.93 


12.70 


-4.65 


BK 


6187.9870 


3.9430 


-1.720 


8.00 


-4.88 


12.72 


-4.57 


FMW 


6191.5580 


2.4330 


-1.417 


81.27 


-5.05 


86.37 


-4.77 


BWL 


6219.2790 


2.1980 


-2.433 


28.82 


-4.89 


37.39 


-4.64 


FMW 


6230.7220 


2.5590 


-1.281 


83.77 


-4.93 


96.11 


-4.66 


FMW 


6246.3170 


3.6020 


-0.733 


58.95 


-5.00 


68.34 


-4.79 


BKK 


6252.5540 


2.4040 


-1.687 


60.75 


-4.97 


76.18 


-4.67 


FMW 


6265.1310 


2.1760 


-2.550 






33.28 


-4.62 


FMW 


6335.3280 


2.1980 


-2.177 


35.89 


-5.02 


46.94 


-4.76 


BWL 


6336.8230 


3.6860 


-0.856 


50.37 


-4.93 


60.57 


-4.72 


BK 


6411.6470 


3.6540 


-0.595 


65.13 


-5.01 


78.87 


-4.62 


BKK 


6419.9420 


4.7330 


-0.240 


44.53 


-4.83 


56.25 


-4.61 


FMW 


6421.3490 


2.2790 


-2.027 


54.10 


-4.83 


63.30 


-4.61 


FMW 


6430.8440 


2.1760 


-2.006 


56.24 


-4.91 


62.60 


-4.73 


FMW 


6496.4660 


4.7950 


-0.570 






33.86 


-4.56 


FMW 
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Table 4: continued. 



Element 






HD 37594 


HD32115 




Wavelength 


^excit 


logg/ 


EQW 


abundance 


EQW 


abundance 


Ref logg/ 


A 


eV 




mA 


dex 


mA 


dex 




6677.9850 


2.6920 


-1.418 






79.98 


-4.67 


BWL 


6750.1500 


2.4240 


-2.621 


18.20 


-4.78 


20.09 


-4.64 


FMW 


6810.2570 


4.6070 


-0.986 


15.00 


-4.79 


18.84 


-4.61 


BWL 


6858.1450 


4.6070 


-0.930 


12.93 


-4.92 


21.44 


-4.60 


BWL 




Average 




-4.91+0.09 (75) 


-4.63+0.07 (88) 




Fen 
















4178.8620 


2.5830 


-2.500 


S 138.9 


-4.80 


S 146.7 


-4.53 


T83av 


4416.8300 


2.7780 


-2.410 


126.64 


-4.93 


136.59 


-4.67 


BSScor 


4491.4050 


2.8560 


-2.700 


116.13 


-4.75 


127.01 


-4.51 


KK 


4508.2880 


2.8560 


-2.250 


141.01 


-4.70 


148.78 


-4.47 


T83av 


4515.3390 


2.8440 


-2.450 


129.10 


-4.75 


144.16 


-4.51 


T83av 


4520.2240 


2.8070 


-2.600 


130.18 


-4.60 


139.43 


-4.41 


T83av 


4541.5240 


2.8560 


-2.790 


94.57 


-4.98 


108.72 


-4.71 


BSScor 


4576.3400 


2.8440 


-2.920 


92.81 


-4.85 


106.02 


-4.56 


T83av 


4620.5210 


2.8280 


-3.240 


63.16 


-4.93 


79.48 


-4.65 


T83av 


4731.4530 


2.8910 


-3.000 


S 83.0 


-4.91 


94.20 


-4.69 


BGHR 


4993.3430 


2.8070 


-3.640 


S 50.0 


-4.83 


S 52.0 


-4.66 


T83av 


5132.6690 


2.8070 


-3.980 


20.63 


-4.98 


29.97 


-4.69 


BSScor 


5197.5770 


3.2300 


-2.100 


133.73 


-4.80 


S 141.0 


-4.58 


KK 


5234.6250 


3.2210 


-2.230 


133.14 


-4.65 


143.70 


-4.43 


HLGN 


5264.8120 


3.2300 


-3.120 


56.62 


-4.85 


S74.9 


-4.62 


T83av 


5284.1090 


2.8910 


-2.990 






94.64 


-4.71 


BSScor 


5325.5530 


3.2210 


-3.120 


S 66.6 


-4.82 


S 82.6 


-4.52 


BSScor 


5414.0730 


3.2210 


-3.540 


S24.3 


-4.96 


S48.6 


-4.55 


T83av 


5425.2570 


3.1990 


-3.160 






68.70 


-4.65 


BSScor 


5525.1250 


3.2670 


-3.950 






S 16.4 


-4.68 


HLGN 


5534.8470 


3.2450 


-2.730 


82.53 


-4.88 


101.58 


-4.56 


BSScor 


5591.3680 


3.2670 


-4.590 


S2.8 


-4.98 


5.71 


-4.55 


RU 


5627.4970 


3.3870 


-4.130 


S 13.0 


-4.70 


S 15.0 


-4.54 


T83av 


6084.1110 


3.1990 


-3.780 


21.80 


-4.85 


31.07 


-4.59 


BSScor 


6113.3220 


3.2210 


-4.110 


S 13.7 


-4.79 


S 19.0 


-4.57 


BSScor 


6149.2580 


3.8890 


-2.720 


51.68 


-4.82 


S67.0 


-4.57 


BSScor 


6247.5570 


3.8920 


-2.310 


89.16 


-4.82 


106.52 


-4.57 


BSScor 


6369.4620 


2.8910 


-4.160 


20.58 


-4.80 


25.30 


-4.63 


BSScor 


6383.7220 


5.5530 


-2.210 


13.74 


-4.80 


19.23 


-4.62 


M 


6416.9190 


3.8920 


-2.790 


50.76 


-4.85 


S 63.6 


-4.60 


M 


6432.6800 


2.8910 


-3.520 


S48.0 


-4.88 






T83av 


6456.3830 


3.9030 


-2.100 






S 120.0 


-4.49 


BSScor 


7224.4870 


3.8890 


-3.240 


25.60 


-4.72 


37.85 


-4.55 


T83av 


7449.3350 


3.8890 


-3.090 


S21.6 


-5.05 


33.22 


-4.74 


HLGN 


7479.6930 


3.8920 


-3.680 


8.92 


-4.85 


14.04 


-4.60 


BSScor 


7515.8310 


3.9030 


-3.460 






26.44 


-4.54 


T83av 


7711.7230 


3.9030 


-2.500 


62.34 


-4.90 






T83av 




Average 




-4.84+0.10(32) 


-4.59±0.08 (35) 
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